Half-Heusler semiconductors as piezoelectrics 
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We use a first-principles rational-design approach to demonstrate the potential of semiconducting 
half-Heusler compounds as a previously-unrecognized class of piezoelectric materials. We perform 
a high-throughput scan of a large number of compounds, testing for insulating character and cal- 
culating structural, dielectric, and piezoelectric properties. Our results provide guidance for the 
experimental realization and characterization of high-performance materials in this class that may 
be suitable for practical applications. 
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One of the central challenges in materials science is 
the design of multifunctional materials, in which large 
responses are produced by applied fields and stresses. A 
rapidly developing paradigm for the rational design of 
such materials is based on the first-principles study of 
a large family of materials. First principles calculations 
of structure and properties are used to explore the mi- 
croscopic origins of the functional properties of interest, 
and this information is used to guide the computational 
screening of many compounds, in both equilibrium and 
metastable structures, to identify promising candidate 
systems. A prototypical example is the optimization of 
piezoelectricity in perovskite oxides. The piezoelectric 
coefficients, describing the strain induced by an applied 
electric field or equivalently a voltage difference induced 
by applied stress, can be readily computed from first prin- 
ciples [TJ |5] • Screening of a large number of systems has 
identified candidate materials for high-performance actu- 
ator and sensor applications [3l 0] . 

Another large family of materials is that of the ABC 
half-Heusler compounds (MgAgAs structure type, also 
called semi-Heusler or Juza-Nowotny compounds [S]), 
with almost 150 distinct compounds reported in the Inor- 
ganic Crystal Structural Database (ICSD) [5J. The half- 
Heusler structure, shown in Fig. n] has FAZm symmetry 
and can be viewed as a rocksalt lattice formed from two 
of the three constituent atoms (at Wyckoff positions 4a 
and 4b) , with the third atom filling half of the tetrahe- 
dral interstitial sites (either Wyckoff position 4c or 4d); 
it is related to the A^BC Heusler structure by removal of 
one of the A sublattices, and can alternatively be viewed 
as a stuffed AC or BC zincblende structure [7]. 

Following an initial wave of interest stimulated by the 
observation of half-metallic ferromagnetic behavior in 
half-Heusler compounds [S] , there has been a resurgence 
of interest in these compounds as materials that can dis- 
play topological properties [HI HO] or be tailored for uses 
as diverse as components in spintronic devices [TT] and 
high-performance thermoelectrics p~2l414j . The semicon- 
ducting half-Heusler compounds are of particular inter- 
est [T51 US]- As insulators, these can exhibit functional 



FIG. 1. (Color online) The ABC half-Heusler structure: A 
(green) and B (orange) are arranged in a rocksalt lattice, with 
the tetrahedral coordination of C (blue) by A shown. 



properties associated with electric polarization, but these 
properties have received very little attention. In fact, 
while the FAZm space group of the half-Heusler struc- 
ture allows a nonzero piezoelectric response, no measure- 
ments of piezoelectricity in these systems have yet been 
reported in the literature. 

In this paper, we use first-principles methods to pre- 
dict the piezoelectric response and related properties of 
half-Heusler compounds. We present these predictions 
first for compounds already reported in the half-Heusler 
structure, and then perform a high-throughput analysis 
of a much larger set of candidate combinations, identi- 
fying high-performance compounds for practical appli- 
cation. We find that half-Heusler compounds exhibit a 
wide range of shear piezoelectric constants; the highest 
values (found for as-yet hypothetical compounds) are well 
above di4 = 200pC/N. This compares well to known 
piezoelectric oxides such as PZT with d^ m 300pC/N 
and ZnO with d^ sa lOpC/N. In addition, the diver- 
sity of possible combinations could allow for other de- 
sirable functional properties to couple to piezoelectricity. 
Through targeted synthesis, which might include compo- 
sitional substitution, epitaxial growth or artificial struc- 
turing, half-Heusler compounds could thus be developed 
as a valuable class of piezoelectric materials, much like 
the perovskite oxides. 



First-principles calculations are performed with the 
ABINIT package [17HT9] using the local density approxi- 
mation (LDA) and an 8x8x8 Monkhorst-Pack sampling 
of the fee Brillouin zone [5U]. Optimized [3T] designed 
non-local [53] norm-conserving pseudopotcntials, gener- 
ated using the OPIUM code [53], are employed, with a 
plane- wave cutoff of 25 Ha. 

For the high-throughput search, we consider combi- 
nations of three distinct elements ABC. We limit the 
search to combinations with a total of 8 s and p valence 
electrons, since we expect this to improve the likelihood 
of band gap formation [7] . We also require that there be 
at least one and at most two p-block elements among the 
three constituents, with the remainder coming from the s 
and/or d blocks of the periodic table. We use Roman nu- 
merals to denote the valences of the constituent elements 
(including d electrons) to classify the combinations into 
families with 8 valence electrons (I-I-VI, I-II-V, I-III-IV, 
II-II-IV, and II-III-III), 18 valence electrons (XI-I-VI, XI- 
II-V, XI-III-IV, I-XII-V, II-XII-IV, III-XII-III, X-II-VI, 
X-III-V, and X-IV-IV), or 28 valence electrons (X-XII- 
VI, XI-XI-VI, XI-XII-V and XII-XII-IV). Among these 
families, we consider members constructed from the se- 
lections I=(Li, Na, K), II=(Be, Mg, Ca, Sr, Ba), III=(B, 
Al, Ga, In, Sc, Y), IV=(C, Si, Ge, Sn, Pb, Ti, Zr, Hf), 
V=(N, P, As, Sb, Bi), VI=(0, S, Se, Te), X=(Ni, Pd, Pt), 
XI=(Cu, Ag, Au), and XII=(Zn, Cd). This generates a 
total of 987 candidate combinations to be searched. 

We first consider the 38 combinations in our search 
set that have been experimentally reported in the ICSD 
in the half-Heusler structure. For each combination, we 
optimize the lattice constant for each of the three struc- 
tural variants ABC, ABC, and ABC, where the under- 
score indicates the unique element that is tetrahedrally 
coordinated by the other two elements. First-principles 
results show that 27 of them, listed in Table [ij are in- 
sulating. The predicted lowest-energy structural variant 
is indicated in the first column. We find that for the 
five compounds (LiZnAs, AuScSn, NiScSb, PdScSb and 
PtYSb) for which a refined structure, including R value 
and temperature factors, is available in ICSD, our pre- 
diction agrees with the experimentally observed variant. 
Most of the computed equilibrium lattice constants given 
in Table IT] are in excellent agreement with experiment. 
There are a few exceptions, LilnSi, LiGaSi and LiZnP, 
for which the DFT values, confirmed by independent all- 
electron calculations, differ significantly from the exper- 
imentally measured values; further experimental investi- 
gation of these cases is warranted. The computed gaps 
range from 0.07 eV for AuYPb and LilnSi to 1.55 eV for 
LiMgAs; measured values are expected to be higher given 
the well-known tendency of DFT to underestimate band 
gaps. 

For each compound, we perform a linear-response 
calculation using density-functional perturbation theory 
(DFPT) [2J], as implemented in ABINIT, to compute the 



electronic dielectric constant e^, dynamical charges and 
zone-center phonon frequencies and eigenvectors, from 
which we obtain the zero-stress static dielectric constant 
eo [IS] , reported in Table [I] Calculations of the strain re- 
sponse [TJ[5B] yield the C44 elastic constant and eu piezo- 
electric coefficient, also reported in Table [TJ These re- 
sults determine du — 614/644 and the electromechanical 
coupling coefficient ku, the conventional figure of merit 
for piezoelectric performance, as ku = \eiA/^/C~44^e^, 
where £f s is the permittivity of free space 1 . We note 
that C44, especially when it is very small, is difficult to 
calculate precisely, the more so as it is also sensitive to 
the value of the lattice parameter. There is a correspond- 
ing uncertainty in du. However, ku is less affected as 
it depends on C44 only as I/1/C44. Also, we note that 
eoo and therefore eo is generally overestimated in DFT, 

— 1/2 

and fci4 oc e , our calculated values of ku are likely 
to be underestimates of the electromechanical coupling 
coefficients that could be achieved experimentally. 

The computed values of eo, £44, eu, du and ku are 
reported in Table [IJ The values of the piezoelectric con- 
stant ei4 range from 0.01 to 0.81 C/m 2 . Seventeen com- 
pounds have ei4 > 0.16 C/m 2 , the experimentally mea- 
sured piezoelectric coefficient of GaAs. The highest pre- 
dicted values are for NiTiSn, PtTiSn and LiGaSi. These 
three compounds also have the highest values of elec- 
tromechanical coupling ku, arising from their large ei4. 
LiMgP, LiZnP, and LiZnAs have similarly high fc 14 val- 
ues, with lower values of ei4 compensated by lower values 
of e . 

It is remarkable that no piezoelectric response data for 
any half-Heusler compound has yet been reported. With 
a single-crystal sample of sufficiently low conductivity, 
the piezoelectric coefficient should be readily measurable 
for most if not all of these compounds. Moreover, mea- 
surements of the dielectric response and elastic coeffi- 
cients, which also have not been reported to date, would 
provide an additional test of these theoretical predictions 
and a more complete characterization of the polarization- 
related properties of these otherwise much-studied com- 
pounds. 

Next, we consider the properties of the full set of 987 
hypothetical and real ABC combinations identified ear- 
lier for study. As before, for each combination we opti- 
mize the lattice constant for each of the variants ABC, 
ABC, and ABQ. Choosing the variant having the lowest 
total energy, we determine whether our LDA calculations 
predict it to be insulating. Of the 987 combinations, we 
find 371 insulators having either 8 or 18 valence electrons, 
while the compounds containing 28 valence electrons are 
all found to be metallic. 

For the insulators, we perform linear-response calcula- 
tions using ABINIT as described above, except that we 
also compute the C\\ and Cyi elastic constants. This 
allows us to screen for local elastic stability by requiring 
that Cn+2Ci2 > 0, C44 > and C n -C 12 > 27 . Fur- 
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LiMgP 


6.02 


5.95 


1.51 


0.32 


7.6 


0.42 


0.15 


12 




LiMgAs 


6.19 


6.17 


1.55 


0.26 


6.9 


0.37 


0.13 


12 




LiMgBi 


6.74 


6.71 


0.62 


0.07 


2.8 


0.26 


0.04 


16 




LiZnN 


4.87 


4.77 


1.21 


0.11 


1.3 


0.84 


0.03 


13 
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LiZnP 


5.78 


5.59 


1.15 


0.44 


6.6 


0.67 


0.15 


14 
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LiZnAs 


5.94 


5.87 


1.12 


0.43 


8.1 


0.54 


0.16 


15 




LiCdP 


6.09 


5.95 


0.76 


0.14 


3.2 


0.46 


0.05 


17 
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AuScSn 


6.42 


6.39 


0.19 


0.17 


2.7 


0.62 


0.05 


20 




AuYPb 


6.73 


6.65 


0.07 


0.44 


9.1 


0.48 


0.13 


27 




LiGaSi 


6.09 


5.81 


0.08 


0.56 


8.6 


0.65 


0.15 


22 




LilnSi 


6.68 


6.18 


0.07 


0.08 


1.6 


0.51 


0.03 


23 




NiTiSn 


5.94 


5.89 


0.44 


0.81 


12.9 


0.62 


0.19 


30 




NiZrSn 


6.11 


6.08 


0.52 


0.19 


2.8 


0.70 


0.05 


26 




NiHfSn 


6.08 


6.00 


0.41 


0.28 


3.4 


0.80 


0.07 


24 




PdZrSn 


6.32 


6.27 


0.49 


0.01 


0.1 


0.68 


0.01 


24 




PdHfSn 


6.30 


6.20 


0.38 


0.06 


0.7 


0.75 


0.01 


22 




PtTiSn 


6.16 


6.15 


0.82 


0.64 


9.5 


0.67 


0.16 


25 




PtZrSn 


6.32 


6.32 


1.06 


0.19 


2.5 


0.77 


0.05 


20 




PtHfSn 


6.31 


6.25 


0.99 


0.23 


2.7 


0.83 


0.06 


18 




NiScSb 


6.06 


6.06 


0.29 


0.15 


2.3 


0.66 


0.04 


21 
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NiScBi 


6.19 


6.15 


0.22 


0.20 


3.6 


0.57 


0.06 


24 




NiYSb 


6.31 


6.28 


0.34 


0.20 


3.6 


0.55 


0.06 


20 




NiYBi 


6.41 


6.36 


0.26 


0.14 


2.9 


0.49 


0.05 


22 




PdScSb 


6.31 


6.27 


0.28 


0.03 


0.4 


0.59 


0.01 


20 




PdYSb 


6.53 


6.47 


0.26 


0.23 


4.2 


0.53 


0.08 


18 




PtScSb 


6.31 


6.31 


0.71 


0.07 


1.0 


0.66 


0.02 


19 




PtYSb 


6.54 


6.51 


0.58 


0.19 


3.3 


0.58 


0.06 


19 
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TABLE I. Properties of experimentally synthesized insulat- 
ing half-Heuslers, grouped by family. Experimental lattice 
constants (in A) are from the ICSD. Also presented are the 
theoretical lattice constant a, band gap -Egap (eV), piezoelec- 
tric coefficients ei4 (C/m 2 ) and di4 (pC/N), elastic constant 
C44 (10 11 Pa), electromechanical coupling coefficient fci4, and 
zero-stress static dielectric constant eo. 



FIG. 2. (Color online) a) Cell parameters (in A) of insulating 
ABC combinations as well as their range of band gaps (in 
eV) are depicted as open black circles, b) Electromechani- 
cal coupling factor k-n of insulating ABC combinations as a 
function of piezoelectric constant ei4. Known combinations, 
from Table III are highlighted as filled red circles. 



thermore, we calculate phonon frequencies at three addi- 
tional high-symmetry points (X, L and W) and eliminate 
combinations which exhibit any unstable modes. This re- 
duces the combinations further from 371 down to 312. 

We thus arrive at 312 combinations that are predicted 
to be insulating and locally stable in the lowest-energy 
variant of the three possible ABC half-Heusler struc- 
tures. Since DFT tends to underestimate band gaps, 
we expect that the actual fraction of insulating struc- 
tures will be slightly higher than our calculations would 
indicate. The computed band gaps and lattice param- 
eters for these 312 compounds are shown in Fig. [2ja). 
We expected no particular correlation between lattice pa- 
rameter and band gap, and indeed we find none. Both 
quantities are rather broadly distributed, suggesting that 
there could be considerable flexibility in choosing mate- 
rials over a substantial range of desired gap or lattice 
constant. 

For these 312 compounds, we compute e^, C44, en, 
du and fci4, using the same methods as before. To give 
a sense of how the range of properties in the full set of 



known and hypothetical compounds compares with that 
of the subset of known compounds, we present a scatter 
plot of fc 14 vs log e 14 in Fig. [2jb). It can be seen that 
there are hypothetical compounds with kn values well 
above those of known compounds (some have fci4 close 
to one). The twenty-seven compounds with the highest 
values of fci4 are listed in Table [TTJ 

As a guide for experimental investigation of piezoelec- 
tric half-Heusler compounds, we highlight a selection of 
candidate compounds chosen according to relevant prac- 
tical considerations. We filter the list to eliminate com- 
pounds with toxic (Pb, Cd, As) or expensive (Be, Pd, 
Pt, Au, In) elements. In addition, we require a DFT 
band gap above 0.4 eV to favor low sample conductiv- 
ity. Lastly, to promote cation ordering into the lowest- 
energy variant, we require AE > 0.15 eV, where AE is 
the difference in energy between the lowest-energy and 
next-lowest-energy variants. Those entries in Table II 
that satisfy these criteria are shown in boldface. Of 
the 100 compounds that satisfy these criteria, 79 have 
ei4 > 0.16 C/m 2 , the experimentally measured value of 
GaAs. 



ABC 


a 


-C/gap 


ei4 


di4 


C44 


fcl4 


eo 


KMgP 


6.71 


1.04 


1.19 


346.28 


0.03 


0.89 


13 


LiNaS 


6.10 


3.12 


0.92 


208.22 


0.04 


0.86 


8 


LiSrAs 


6.84 


1.93 


0.79 


286.27 


0.03 


0.85 


10 


LiNaSe 


6.39 


2.70 


0.75 


227.09 


0.03 


0.84 


8 


NaMgN 


5.43 


0.72 


1.50 


165.23 


0.09 


0.83 


13 


NaBeP 


5.59 


2.06 


1.02 


192.98 


0.05 


0.82 


11 


KMgAs 


6.92 


0.67 


1.06 


205.75 


0.05 


0.81 


13 


KScPb 


7.35 


0.42 


1.21 


330.83 


0.04 


0.80 


25 


MgBaSi 


7.18 


0.37 


1.13 


290.51 


0.04 


0.80 


21 


BeBaSn 


6.86 


0.16 


0.49 


1036.54 


0.01 


0.80 


33 


KScSn 


7.31 


0.42 


1.19 


294.69 


0.04 


0.79 


24 


KBeSb 


6.51 


0.81 


0.72 


266.50 


0.03 


0.79 


13 


MgBaGe 


7.21 


0.06 


1.10 


266.49 


0.04 


0.78 


21 


AgSrAs 


6.72 


0.34 


1.10 


163.79 


0.07 


0.77 


14 


NaZnP 


6.02 


0.76 


1.67 


163.79 


0.14 


0.77 


16 


NaBSi 


5.43 


0.44 


0.94 


121.79 


0.06 


0.73 


13 


KBaN 


6.48 


0.69 


0.91 


174.57 


0.05 


0.70 


18 


NaZnAs 


6.26 


0.16 


1.44 


98.70 


0.15 


0.69 


18 


KYGe 


7.13 


0.16 


1.06 


178.85 


0.06 


0.68 


25 


CaBAl 


5.81 


0.37 


1.61 


99.88 


0.16 


0.68 


21 


KMgSb 


7.30 


0.69 


0.87 


108.05 


0.08 


0.66 


14 


KYSi 


7.11 


0.14 


1.03 


167.49 


0.06 


0.66 


25 


CuSrAs 


6.39 


0.36 


0.69 


179.58 


0.04 


0.65 


19 


NaKO 


5.90 


1.88 


0.60 


80.78 


0.07 


0.64 


8 


KBeBi 


6.62 


0.34 


0.69 


127.24 


0.05 


0.62 


16 


KAISi 


6.57 


0.54 


0.98 


101.96 


0.10 


0.62 


18 


ZnYB 


5.78 


0.33 


1.84 


95.35 


0.19 


0.62 


32 



TABLE II. Top half-Heusler candidate compounds, ranked 
according to predicted electromechanical coupling coefficient 
ki4. Those that remain after applying the filtering criteria, 
as described in the text, are in boldface. Units as in Table 1. 
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LiAgTe 6.38 1.02 0.80 28.60 0.28 0.36 17 



LiCuTe 


6.06 


1.57 


0.38 


8.14 


0.38 


0.16 


14 


LiScC 


5.23 


0.90 


0.34 


7.73 


0.44 


0.12 


19 


LiCuS 


5.43 


1.32 


0.22 


4.04 


0.55 


0.10 


10 


LiMgSb 


6.59 


1.30 


0.15 


4.94 


0.30 


0.08 


14 


LilnSi 


6.18 


0.07 


0.08 


1.64 


0.51 


0.03 


23 



TABLE III. Piezoelectric properties of the previously unchar- 
acterized combinations from Ref. [28] in which the ground 
state is predicted to be the half-Heusler structure. Units as 
in Table I. 



Further investigation of the hypothetical half-Heusler 
piezoelectrics hinges on the possibility of experimentally 
realizing the desired compounds in the half-Heusler struc- 
ture. Additional information about bulk equilibrium 
ABC phases can be obtained from the ICSD. In the cases 
where no phase is reported, it could be that either no 
stable bulk phase exists with that stoichiometry, or sim- 
ply that the relevant composition has not been studied. 
ABC combinations are also reported in several struc- 
tures other than the half-Heusler structure. Specifically, 



we find that twelve combinations listed in Table |TT] are 
reported in ICSD with other structures. Nine (KMgP, 
LiNaS, MgBaSi, MgBeGe, NaZnP, NaZnAs, KMgSb, 
NaKO, and KMgAs) are reported in the P4/nmm PbCIF 
(also called Cu2Sb) structure, one (LiNaSe) in the Pnma 
MgSrSi structure, and two (AgSrAs and CuSrAs) in the 
PQ^/mmc ZrBeSi structure. One (NaZnAs) is reported 
in both PA/nmmm and FmZm structures. 

First-principles calculations of the total energy of alter- 
native structures can also be used both to predict ground 
state structures and to identify systems with piezoelec- 
tric structures as low-energy alternatives that would be 
suitable candidates for stabilization through composi- 
tional substitution or epitaxial strain. In a recent first- 
principles study [28] of the stability of ABC compounds, 
sixteen new combinations are predicted to form in the 
half-Heusler structure. Of these, fifteen are in our search 
set, and we find that six are insulating. Our results for 



these compounds are given in Table III where the com- 
pounds are ranked by feu. It should be noted, however, 
that 462 of the 987 combinations in our search set were 
not considered in the stability study; if they had been, 
the number of compounds in Table |HI| would likely be 
larger and span a wider range of property values. In any 
case, absence from this list does not preclude the synthe- 
sis of the half-Heusler structure as a metastable phase 
with appropriate processing. For example, our calcula- 
tions on KBaN indicate that the half-Heusler structure is 
the second-lowest in energy among seven structures stud- 
ied, and is only 0.11 eV higher than the lowest-energy 
(Pnma) structure; such a case could be promising can- 
didate for directed synthesis as a metastable phase. 

In summary, we have used a first-principles rational- 
design approach to demonstrate semiconducting half- 
Heusler compounds as a previously unrecognized class 
of piezoelectric materials. We have presented these pre- 
dictions first for compounds already reported in the half- 
Heusler structure, and then for a much larger set of can- 
didate combinations that were generated and screened 
for high performance via a high-throughput analysis. 
We hope that our results will provide guidance for 
the experimental realization and further investigation of 
high-performance materials suitable for practical appli- 
cations. We also suggest that the combination of piezo- 
electric properties with other characteristic properties of 
Heuslers, especially magnetic properties, may offer invit- 
ing avenues for further development of multifunctional 
materials. 
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